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We analyze two types of full-circle angle 
calibrations: a simple closure in which a 
single set of unknown angular segments is 
sequentially compared with an unknown 
reference angle, and a dual closure in 
which two divided circles are simulta- 
neously calibrated by intercomparison. In 
each case, the constraint of circle closure 
provides auxiliary information that (1) en- 
ables a complete calibration process without 
reference to separately calibrated reference 
artifacts, and (2) serves to reduce measure- 
ment uncertainty. We derive closed-form 
expressions for the combined standard un- 
certainties of angle calibrations, following 
guidelines published by the International 
Organization for Standardization (ISO) and 



NIST. The analysis includes methods for 
the quantitative evaluation of the standard 
uncertainty of small angle measurement us- 
ing electronic autocollimators, including the 
effects of calibration uncertainty and air 
turbulence. 
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1. Introduction 



Full-circle calibrations of the angular divisions of 
polygon mirrors, indexing and rotary tables, and angle 
encoders can be accomplished without reference to sep- 
arately calibrated reference artifacts using comparator 
techniques together with the principle of circle closure. 
The latter is a natural conservation law for plane angle, 
known since the time of Euclid, expressing the fact that 
the sum of the angles around any point in a plane equals 
27Trad (360°). Should a circle be divided into n angular 
segments A 1,^2, . . . , A„, and the difference between 
each segment and an unknown reference angle X be 
measured, closure provides a constraint on the data that 
enables a complete solution for all n + I unknowns. 
Circle closure is one of a number of self-proving com- 
parator techniques employing multiple measurements 
together with suitable rearrangements of the compo- 
nents of a measurement system. A review of such tech- 
niques and their applications in dimensional metrology 
is given in Ref. [1]. 



In this paper we analyze two types of circle closure: 
a simple closure, just described, and a dual closure 
wherein two artifacts, such as a pair of indexing tables 
or an indexing table and an optical polygon, are inter- 
compared in such a way that each is calibrated in the 
process. The conceptual ideas here are well known and 
widely applied [2-9]; our major goal is to derive explicit 
expressions for the resulting measurement uncertainties 
in a manner consistent with guidelines published by the 
International Organization for Standardization (ISO) and 
NIST [10-11]. In the process, we show how closure 
constraints serve to reduce uncertainty in an interesting 
way as a result of a basic observation: the result of a 
measurement provides information not only about the 
measured quantity (i.e., the defined quantity subject to 
measurement, or the measurand), but also about any 
other quantities with which the measurand shares a 
functional relationship. 
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2. Simple Closure 

In a simple closure calibration, the n angles of interest 
are independently compared to an unknown reference 
angle X using a difference measurement procedure. We 
consider the common case where the angles of interest 
Ai, A2, . . . , An and the comparison angle X are nomi- 
nally equal to Itt/h radians. Denoting the deviations 
from the nominal value hy a\, a2, . . . , a„ and x, respec- 
tively, we have: 



Ai = lirln + Ui, i = 1, 2, . . . , n, 
X = Itt In + X, 



(1) 



where the deviations are assumed to be small angles that 
can be conveniently measured using an accurate 
autocollimator. With these definitions we have A, - Z = 
a, - X , / = 1 , 2, . . . , n and the constraint of circle closure 
becomes: 



where nii = 62 - 6i. We then re-position the dihedral 
mirror to bring the first face back to a position near the 
autocollimator null and repeat the measurement se- 
quence. Continuing in this fashion, we sequentially 
compare the angle X with each of the indexing table 
intervals, yielding the complete set of measurement 
equations: 



fll - X = OTi 

a2 - X = m2 
a^ - X = m^ 



(4) 



Here we have n linear equations in the n + 1 unknown 
quantities {ci, a2, . . . , a^; x] , constrained by the closure 
relation, Eq. (2). Adding Eqs. (4) yields: 



2_i cii - nx = 2j rrii , 



(5) 



2t7 =ZjAi=2j i^Tr/n + a,) 



1=1 1=1 



= n ■ (2TT/n) + 2j <^i = 2t7 + 2j ^^ 

i=l 1=1 

or, finally: 

n 

E«. = 0. (2) 



The physical realization of the reference angle X will 
depend on the chosen measurement procedure. Two 
particular procedures that lead to identical sets of mea- 
surement equations are (a) the calibration of an n -sided 
optical polygon using two autocollimators, and (b) the 
calibration of an indexing table, at a subset n of its 
discrete positions, using a dihedral mirror or a segment 
of an optical polygon. In the first case, X is the angle 
between the optical axes of the autocollimators; for the 
indexing table calibration, X is the angle between the 
two mirror normals. 

Figure 1 shows the setup for a simple closure calibra- 
tion of an n -position indexing table using a dihedral 
comparator mirror (Fig. la) and an electronic autocolli- 
mator. The mirror is aligned so as to bring one of its 
faces near the null of the autocollimator, and a reading 
is taken (Fig. lb). The table is then indexed through 
angle Ai, rotating the second mirror face into near- 
normal incidence, and a second autocollimator reading 
taken (Fig. Ic). Denoting the measured angles by 9i and 
62, we then have, for a suitable sign convention: 



and since the first sum vanishes we have the comparator 
angle deviation x , as an explicit function of the measure- 
ment data: 



1 

X = --^ rrii . 



(6) 



Now using Eq. (6) for x and the set of measurement 
relations Eq. (4), we have for the deviation of the A:th 
segment of the indexing table: 

aii = nik + X 

1 " 

n f. 



rrik- 



nik 1 



2 >ni 



n- 1 



n J n 



- 2 '"; 



(7) 



02 - d\ 



(3) 



We note from the results, Eqs. (6) and (7), that the 
angular units for closure-based angle metrology are 
fixed by the scale of the small-angle measuring system, 
assumed here to be a calibrated autocollimator. The set 
of Eqs. (4) is invariant under the change of scale 
{a,, X, m,- } — > {aa,, ax, arrii }, / = 1, 2, . . . , n where a 
is an arbitrary scale factor. Thus, while closure obviates 
the need for separately calibrated comparator artifacts, 
the achievable uncertainty is limited by the degree to 
which the SI derived unit of plane angle (the radian) is 
realized by the small-angle measuring system. 
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Fig. 1. Calibration of an n -position indexing table using a dihedral compara- 
tor mirror and an autocoUimator. The measurement sequence shown yields the 
first of n difference equations. 



3. Measurement Uncertainty 

In order to evaluate the measurement uncertainty, we 
need the following result from the analysis of the propa- 
gation of uncertainty [10-11]. Suppose that a measur- 
and of interest Y is related to n independent and uncorre- 
lated "influence quantities" zt, k= 1,2, . . . ,n, via the 
functional relation 



u^izk)- The positive square roots of these variances, 
u{zk), k = I, 2, . . . , n, are called the standard uncer- 
tainties of the quantities zu- A first-order Taylor series 
approximation of Y = f{zj, Zi, ■ ■ ■ , Zn) then yields the 
estimated variance of Y: 



"■« = l(l 



u\zk). 



(9) 



Y=f{zi, Z2, ■ ■ ■ , z„) 



(8) 



In a procedure to determine Y, the quantities Zk may be 
directly measured or simply estimated; in either case our 
state of knowledge of these quantities is described by a 
set of probability distributions with associated variances 



with the combined standard uncertainty of Y given by 
u^iY). The quantities df/dzk are known as sensitivity 
coefficients whose squares are essentially weighting 
factors for the individual variances u^{zk)- 
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We emphasize here that the use of probabihty distri- 
butions to characterize our state of knowledge should 
not be taken to imply that the quantities Zk are intrinsi- 
cally random in nature. A dimensional measurement in 
engineering metrology is a deterministic exercise in 
which time and cost constraints determine the quantity 
and quality of the input measurement data. The quanti- 
ties Zk are generally a set of unknown constants whose 
values at the time the data were taken are not known 
exactly. Certain influences, such as the effect of air 
turbulence on the reading of an autocollimator, may 
appear at first sight to fluctuate randomly and unpre- 
dictably. When this happens we often choose to average 
over many readings in order to approximate a mean, or 
equilibrium, value for the measured angle. On the other 
hand, the time-dependent spatial distribution of the air 
density between an autocollimator and a target mirror 
could, in principle, be measured very accurately using 
interferometric techniques, and the result used to calcu- 
late, and correct for, the subsequent refraction error in 
the optical signal. That we normally choose not to do so 
is a matter of economics and not a decision forced on us 
because of any randomness inherent in the measurement 
process. In using a probability distribution to describe 
an experimental parameter, what is "distributed" is not 
the parameter itself but rather our knowledge with re- 
spect to the value it had when we performed the mea- 
surement. 

For the comparator angle x , we then use Eq. (9) to 
compute the variance: 



"'^"^=S(^)'"° 



E 



1 I 2 ^0 ^0 



n I n n 



(12) 



so that the combined standard uncertainty of x is: 



u^{x) 



Mo 



(13) 



Equation (13) expresses what we might intuitively 
expect on the basis of repeated measurements in the 
presence of uncorrected measurement "noise." Since 
the unknown angle x has been the subject of n indepen- 
dent and uncorrected difference measurements, the 
uncertainty of the measured value of x has been reduced 
by V n from the uncertainty of a single measurement. 
Then, using Eq. (10) and the assumed variances 
M^(m,) = Mo, ( = 1, 2, . . . , n, we have: 



uliuk) = 2 VK^y u\mi) 



"-'^' ^ + y,f^>"^ 



,^A«' 



" ( dx\^ 



(10) 



= Mo 



n — \y n — \ 

+ T- 



In Appendix A we show that for an accurately 
calibrated autocollimator measuring angular artifacts of 
high quality (in the sense that measured angles are rea- 
sonably small), the variance M^(m,) may be assumed to 
be constant for all measurements nii, m2, ... , m„. We 
suggest two ways in which this constant value can be 
assigned in practice; for typical laboratory environ- 
ments, it is reasonable to assume that: 



Mo 



(14) 



Thus, for the combined standard uncertainty of Uk we 
have: 



Uciau) = 



1 



■ Mo, 



(15) 



u^(mi) 



213^ 

N 



sUR) 



ul ■ 



constant. 



(11) 



where (3 is a calibration parameter on the order of unity 
and Sn (R) is the experimentally measured (Type A) 
variance of the mean R of a set {Ri, R2, ■ ■ ■ , Rn} of N 
autocollimator readings for a fixed target mirror 
position. 

With this result, Eq. (10) becomes: 



with the same result for the remaining angular devia- 
tions. Equation (15) illustrates an interesting feature of 
calibrations that are constrained by subsidiary condi- 
tions such as circle closure. Given a single measure- 
ment of ak, we would expect an uncertainty equal to Mo, 
w hereas the actual result is less than Mo by a factor of 
V(n - l)/n. Each measurement of one of the deviations 
a, ¥= Uk adds a little to our knowledge of Uk because of 
the closure constraint. As the number of unknown 
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deviations n increases, the factor V(n - \)ln 
approaches unity and the reduction in uncertainty de- 
creases, as we might expect. The closure constraint only 
provides a single piece of information whose effect is 
diluted as the number of unknowns becomes large. 

4. Dual Closure 

Indexing tables are commonly calibrated at a set of n 
positions using an intercomparison technique wherein 
the table of interest (the test table) is mounted concentri- 
cally atop a comparison table capable of indexing to the 
same n positions. Typically one table is rotated clock- 
wise (CW) through a chosen angle, followed by a coun- 
ter-clockwise (CCW) rotation of the second table 
through the same angle, so that a plane mirror target 
mounted on the top table is returned to a nominal zero 
position. Small differences in mirror position before and 
after indexing are measures of differences between pairs 
of nominally equal angular segments of the two index- 
ing tables. These small difference angles are normally 
measured using an accurate autocollimator. 

By the nature of the procedure, the angular position- 
ing errors of both indexing tables are simultaneously 
determined by a complete set of measurement data. 
Therefore, in what follows, we no longer distinguish the 
"test table" from the "comparison table" but simply 
refer to them as the top table (T) and the bottom table 
(B), respectively. We denote the n angular segments of 
interest for the top and bottom tables, respectively, by 
r, = 2ir/M + ti and Bi = 2TT/n + bi,i= 1,2, . . . , n , where 
(bi, ti) are small angular errors. Then for all ( we have 
Bi - Ti = bi - ti, and the circle closure constraints on the 
two tables can be written, by the steps leading to Eq. (2): 

n 



^ h = o 



(16) 



segments Ti and Bi initially positioned as shown in the 
top diagram of Fig. 2. The small dots locate the zero 
index marks on each table. In acquiring measurement 
data, the bottom table always rotates CCW while the top 
table rotates CW, viewed from above. 

Now consider the measurement sequence consisting 
of the following moves: (1 ) rotate the bottom table CCW 
through angle Bi and then (2) rotate the top table CW 
through angle Ti. This sequence is illustrated in the 
middle and bottom diagrams of Fig. 2. The net result 
of the two moves is a CCW rotation of the bottom table 
through a nominal angle of Bi~2tt/3 rad (120°) 
and a rotation of the target mirror through angle 
Bi-Ti = bi- ti, for a suitable set of sign conventions. 
Denoting by mj the difference between the initial and 
final autocollimator readings, we have 

bi - ti =mi. (17) 

We now repeat the measurement sequence two 
more times, yielding two more difference equations 
m2 = b2- t2 and m^ = b^ - t^. At the end of the three 
measurements, the entire system is returned to its initial 
configuration. We then rotate the top table through angle 
Ti ~ 120° CW, re-position the target mirror near the 
autocollimator null, and repeat the set of three two-move 
sequences. Proceeding in this manner yields a complete 
set of 3 X 3 = 9 difference measurements in which the 
three angles of the top table are each compared with the 
three angles of the bottom table. (We note here that this 
is only one of many different complete move sequences 
that yield 3X3 = 9 independent differences between 
pairs of table segments. For any given sequence, the 
nature and numbering of the coefficient matrix and data 
vector defined below will differ in detail from that 
derived here but the resulting uncertainty evaluations for 
the table angle deviations will be identical.) 

The 3X3 indexing table difference equations result- 
ing from a complete set of segment comparisons are 
given by 



5. Explicit Example: « = 3 

In order to demonstrate the calculation of measure- 
ment uncertainty for indexing table calibrations using 
dual closure, we will work out the results for two three- 
position tables in enough detail so that all steps in the 
procedure may be appreciated. We will then be in a 
position to generalize to the n X n case. We thus 
consider the setup shown in Fig. 2, with a pair of three- 
position tables. 

We label the three segments of the top table {Ti, Ti, 
T3} and those of the bottom table {Bi, B2, B3}, with 



bi - ti = m\ 
b2- t2 = m2 
bi- t3 = rrii 
bi - 12 = 11X4 
b2-ti = ms 
bi-t,= me 
b[ - 13 = rrij 
b2-t,= mg 
%-t2 = mg , 



(18) 
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Fig. 2. Dual closure calibration of a pair of 3-position indexing tables. The two-move 
measurement sequence shown yields the first of difference equations. 



where mj, . . . , mg are differences between pairs of 
autocollimator readings. This over-determined set of 
nine linear relations among six unknown quantities has 
no unique solution. This can be seen by noting that 
set of equations is invariant under the additive transfor- 
mations bi -^ bi + K,ti —^ti + K,i= 1,2,. . . ,n, where 
K is an arbitrary constant. The equations are also 
not independent; we see, for example, that 
m^ = m4 + m^ - mi. 

In order to proceed, we combine the measurement 
relations, Eqs. (18), with the two constraint Eqs. (16) 
and express the resultant set of linear equations in 
matrix form: 
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(19) 
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Ax = m. 



(20) 



where A is an 1 1 X 6 coefficient matrix, jc is a 6 X 1 
column vector of indexing table angular deviations: 
x = {bib2b^ f 1 ^2 ^3) and /« is an 1 1 X 1 vector of measure- 
ment data plus constraints. This is an over-determined 
set of 1 1 equations in 6 unknowns. Such data sets are in 
general inconsistent because of measurement errors, 
which means that there is no solution vector jc which 
satisfies all eleven equations simultaneously. This is a 
standard problem in linear algebra and can be solved (in 
the sense of closest approximation) by the method of 
least squares. For a detailed development see, for exam- 
ple, Ref. [12]; in this paper we simply state the result. 
The rank of a matrix is equal to the number of linearly 
independent rows (or columns) of the matrix. Given the 
set of linear equations Ax = m , where the m X n matrix 
A has rank n, then the least-squares solution jc = 
(bi 1)2 h t\ h h) is given by: 



JC = (A^ A)-' A^m, 



(21) 



where A ^ is the transpose of the coefficient matrix A . In 
our 3X3 problem, it can be shown that the rank of A 
is equal to 6, as required, and so we proceed to work out 
the product matrix: 



{A^A\'A 
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(22) 

Then, using jc and m above, with Eqs. (21) and (22), and 
carrying out the matrix multiplication, yields the least- 
squares estimates for the angle deviations of the index- 
ing tables: 



J_ 

27 



[7(oti + OT4 + nil) 



■ 2{m2 + m^ + nis + m^ + m^ + mg)], 
- 1 



(23) 



27 



[7(ot2 + OT5 + mg) 



^^ = YJ ^'^^"^^ + m(, + OT9) 



2(oti + m2 + m^ + m-i + nij + mg)], (25) 

fl = 2^ [-7(OTi + OT6 + Wig) 

+ 2(ot2 + OT3 + OT4 + ms + m? + wtg)], (26) 



fc = 2^ [-7(OT2 + OT4 + OT9) 

+ 2(oti + m^ + nis + m^ + nij + mg)], (27) 



^3 = 2^ [-7(OT3 + OTs + OT?) 

+ 2(oti -H OT2 + ^"4 + w^e + w^s + Wig)], (28) 

Equations (23)-(28) display the explicit functional rela- 
tionships between the indexing table angle deviations 
and the measurement data. The combined standard 
uncertainties can then be calculated by taking the appro- 
priate derivatives according to Eq. (9). For the variance 
in the least-squares estimate of fo, we have, using 
Eq. (23) 



«c^(^.)=|:(0«v.) 



(29) 



As before, the quantities dbildnii are the sensitivity 
coefficients and here we see that measurements OTi, OT4, 
and lUj have greater weight, by a factor of 7/2, than that 
of the remaining measurements. This follows naturally 
since these three measurements depend directly on the 
value of b\ [see Eqs. (18)], while the remaining mea- 
surements depend on bi only indirectly, by way of the 
constraint relations. As in the case of simple closure, the 
measurement data here consists of a set of independent 
autocollimator difference readings and we again make 
the reasonable assumption (see Appendix) that each 
measurement has the same uncertainty: u^inii) 
i = I, 2, . . . , 9. Then Eq. (29) becomes 



Mo, 



19 



u^(b.) = — 2 [3X7^ + 6 X2^]ui = — ui , (30) 



27 



81 



with the combined standard uncertainty given by 



2{m\ + m^ + m^ + m^ + nij + OT9)], 



(24) 



Mc(Ol) = —p— Mo 



(31) 
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The uncertainties of the remaining least-squares 
estimates {b2, b^, fj, h, h) are also given by Eq. (31). 

6. Generalization for Two n -Position In- 
dexing Tables 

The procedure described for the 3X3 indexing table 
intercomparison may be generalized to the case of an 
n X n calibration. The set of all difference measure- 



ments bj-tj, i,j = 1, 2, 



, n between the angular 



segments of the two tables will yield a set of n^ equa- 
tions analogous to Eqs. (18). Using the same angular 
segment labeling conventions, the general problem will 
have an angle deviation vector x with 2n components 
(bi Z?2 . . .b„ti t2. ■ ■ t„), a coefficient matrix A of size 
(n^ + 2) X 2n, and a data vector m of length n^ + 2. 
Using a matrix manipulation code such as MATLAB or 
Mathematical it is not difficult to show by induction 
that A has rank 2n and that Eq. (23), for general n, 
becomes: 



1 

3«^ 



[(3n -2)(mi + ot„+i 



+ OT2n+l + . . . + m„(„_i)+i) 
- 2(OT2 + m-i+ . . . + OT„ + OT„ + 2 + OTn + 3 +. . . +OT2n 

+ m2n + 2+- ■ ■+m„2)], (32) 

with analogous expressions for b2, ■ ■ ■ , b^ and 't\, . . . , 
In. The numbering scheme makes it unwieldy to write 
these expressions explicitly, but the general rule is 
simply 

h = ^2 [(3« - 2) 

X (2(all m's with explicit dependence on bu)) 

- 2 (^(remaining m 's))] , (33) 



and 



Ik = :r^ [-(3« - 2) (2(all m's with explicit 



dependence on 4)) + 2 (^(remaining m's))], (34) 



' Certain commercial equipment, instruments, or materials are identi- 
fied in this paper to foster understanding. Such identification does not 
imply recommendation or endorsement by the National Institute of 
Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the 
purpose. 



where A: = 1 , 2, . . . , n. 

The calculation of the combined standard uncertainty 
proceeds in exactly the same way as that leading to 
Eq. (31), under the reasonable assumption (see 
Appendix A) that the combined standard uncertainty of 
each measurement m; is just u^{mi) = Mq. Examining 
Eq. (32), we see that the quantity multiplied by (3n - 2) 
contains n terms, while the quantity multiplied by 2 
contains n{n - 1) terms. (The total number of m -terms 
is just n^, the number of measurements.) We can then 
write, using Eq. (29) 



Mc^^l) 



3n^ 



[nOn-lf + Anin - 1)]mo 



9n 



(9n - 8)mo' 



\n 9n' 
so that the combined standard uncertainty is 



(35) 



Uc{b\) = 



9n 



1/2 

2 ; Mo . 



(36) 



The same result follows for the combined standard un- 
certainty of the remaining angle deviation estimates. 
Equation (36) i£ quite interesting. If we neglect the term 
in n^, then u^b,) = Uo/vn, which is what we might 
expect intuitively, since the unknown angle Bi has been 
independently compared with the n angles Ti, . . . , T^. 
The factor 8/(9n^) reduces the uncertainty slightly as a 
result of the constraints due to circle closure. Since 
there are only two closure relationships for any n , while 
the number of measurements grows like n^, the reduc- 
tion in combined standard uncertainty becomes 
insignificant for large n; for n = 12 for example, a 
common value, the correction is less than 5 %. 

7. Conclusions 

Angle calibrations using closure constraints are self- 
proving and can be accomplished without the need for 
separately calibrated reference artifacts or the services 
of a standards laboratory. The SI unit of angle is 
transferred by a small angle measuring system, typi- 
cally a calibrated electronic autocollimator. The con- 
straints lead to (usually modest) reductions in the com- 
bined standard uncertainties of the measured angles. 
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These uncertainties can be readily calculated using the 
results derived in this paper. In common laboratory en- 
vironments, the primary source of measurement uncer- 
tainty for a well-characterized autocollimator is turbu- 
lence noise on the optical signal. The minimum 
achievable uncertainty is limited by the degree to which 
such noise can be reduced. 



8. Appendix A. Autocollimator 
surement Uncertainty 



Mea- 



An autocollimator is an optical instrument used to 
measure small rotations of a plane mirror target. The 
output (reading) of a typical electronic autocollimator is 
a number proportional to the angle between the normal 
to the target mirror and the instrument optical axis. 
Since there is no direct physical access to the optical 
axis, a single reading with a fixed target conveys no 
useful information. Autocollimators are therefore used 
in either a differential mode wherein the difference be- 
tween two readings measures target mirror rotation, or 
in a nulling mode where changes in target orientation 
are sensed and used for optical system alignment or 
servo control. 

A first-principles analysis of autocollimator measure- 
ment uncertainty is by its nature very complex. The 
measurand as defined is the angle between the normal to 
a perfectly flat mirror and the propagation vector of a 
perfectly plane wavefront emitted by, and subsequently 
collected and focused by, a perfect optical system onto a 
perfectly linear detector. In practice, all of these features 
are realized only approximately. We are saved by the fact 
that real autocollimators, while useless for absolute 
measurements of angular position, can be accurately 
calibrated to measure small rotations or changes in an- 
gular orientation, even with imperfect optics and target 
mirrors that are not really flat. 

The response of an autocollimator is in general non- 
linear, especially for large angles near the limits of the 
measuring range. For calibrations typical of those con- 
sidered in this work, the setup can be arranged so that 
all measurements occur near the center of the instrument 
range so that the response can be assumed to be linear. 
Calibration data for a particular instrument will serve to 
justify this assumption. For a single reading R of such an 
autocollimator, the sensed angle 9 is given by: 



e=pR 



(Al) 



where ;S is a proportionality constant derived from the 
instrument calibration. Because of refraction caused by 
turbulence in the air path between the autocollimator 
output aperture and the target mirror, the instrument 



reading will fluctuate about a mean wlue. It is common 
measurement practice to minimize this air path, employ- 
ing shielding where possible, and to reduce the effects of 
turbulence by averaging over many readings. 

Now suppose that a small rotation angle is measured 
by averaging A^ autocollimator readings at each of two 
mirror positions {9i, initial, and ft, final) and then tak- 
ing the difference of these averages. Denote the two sets 
of readings by Rij and R2k, j, k = 1, 2, . . . , N. Using 
Eq. (Al), we then have for the measured angles: 

- B ^ 



^^ - Tt Zj ^2k - PRi 



and the measured mirror rotation m is just 



(A2) 



-n _^ 



P 



- O2 - 6, - — 2j R2k-T}2j Rij 



= m(;8, R2U R22, ... ,R2n; Rn, Rm, ... , Rm), (A3) 

where the last expression shows the explicit dependence 
of m on the measurement data. We then compute the 
variance of m, using Eq. (9): 



uXm) = (^)' u\l3) 



N I 
k= I \ 



dm Y 2 
I It 

dR2k. 



u\R2k) 



N , 






(A4) 



From Eq. (A3), the sensitivity coefficients are given by: 



dm 1 



= iS R2J.-h l.Rii = ^2-Ru (A5) 



dp N 



N 



and 



dm _ dm _ (3 



(A6) 
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The variance u^ijS) in Eq. (A4) characterizes the stan- 
dard uncertainty associated with the autocollimator cal- 
ibration; for the variances u\R2k) and u\R\j) we use the 
Type A (statistical) evaluation estimates (sample vari- 
ances): 



1 " 
u\R2k) = s\R2k) = j^—^ E (^^2, - Rlf , 

1 " 

\R,,) = s\R,j) = j^-^ 2 (^1. - Rif , (A7) 



With the results Eqs. (A5)-(A7), Eq. (A4) becomes 



u^im) = iR2 - R,f u\l3) + ^5 lARik) + s\R,j)] 



{R2 - R^' u\P) + ^ [s\R2) + s\R>)] , (A8) 



where s^{R2) = s^{R2k)IN , the experimental variance of 
the mean R2 and similarly for s^{R\). For a particular 
environment, the measurement noise induced by turbu- 
lence is normally stable in the sense that the standard 
deviation of the mean of a sample of A^ readings is 
sensibly constant from sample to^ sample. We therefore 
assume that s^{R2) = s^{R\) = Sn{R) = a constant for each 
sample of A^ readings. We then have for the variance of 
the angle measurement m 



Uc(m) = {R2 



■R,fuM3) + ^sUR) 






P' 



(A9) 



measuring time and susceptibility to thermal drift. 

The dependence of u^{m) on m^ means that each 
measurement in a full-circle angle calibration will, in 
general, have a different combined standard uncertainty. 
In practice, the set of individual variances can be ap- 
proximated by a single constant value with negligible 
increase in overall measurement uncertainty. There are 
two ways to effect this approximation. 

We first note that for an accurately calibrated autocol- 
limator and typical high-quality angular artifacts, the 
measurement uncertainty will often be dominated by 
turbulence noise. To see this, we examine Eq. (A9) using 
a representative set of measurement parameters. A 
modern electronic autocollimator is adjusted by the 
manufacturer to set y8 = 1 , and after calibration we could 
expect a worst-case relative standard uncertainty of 
m(;8)/|;8| < 0.1 % = lO'l In the process of calibrating 
an indexing table with a dihedral mirror matched to the 
nominal angles of interest, the maximum angle that we 
might expect to measure would be on the order of 
25 (xrad (1 (xrad = 0.2"). Thus mV{P)/p^ < 6.3 X 
10^ (xrad^. In a typical measurement setup at NIST, 
with an air path of several centimeters between the auto- 
collimator and the target mirror, we average over several 
hundred autocollimator readings per measuring position 
and_commonly measure turbulence noise values near 
5^(^)^=0.25 (xradl 

Comparing these numbers, we see that in such exper- 
imental conditions the measurement uncertainty due to 
turbulence noise is larger than that due to calibration 
error by more than two orders of magnitude. Thus we 
can, with negligible error, set ul{m) ~ {213^/N) Sn{R) = 
Mo, a constant for all measurements, with the actual 
value assigned to Mq determined by locally measured 
turbulence noise characteristics. 

A second, more conservative way of assigning a con- 
stant uncertainty is to use the maximum value computed 
from Eq. (A9), based on the measurement data. Setting 
w^MAx = MAX{mi, 1712, . . . m„ }, we then fix a constant 
value for Uo by setting 



Mo 



UciniMAx) = muAx 



2 u\§) 20^ - 
B^ N "^ ' 



(AlO) 



where we have used Eq. (A3) in the form {R2 - Rif = 
rn^ip^. Equation (A9) is the central result for an autocol- 
limator measurement with a negligible linearity error. 
The first term on the right-hand side arises from calibra- 
tion uncertainty, as measured by the square of the 
relative standard uncertainty M(y8)/|;8|. The second 
term, arising from turbulence noise in the autocollima- 
tor signal, can be reduced by increasing the number of 
autocollimator samples A^ at the cost of increased 



This will result in overestimating the variances and com- 
bined standard uncertainties of some of the measured 
angles, but the turbulence component will still account 
for most of the uncertainty in practice. Either of the 
techniques just described can be used for quantitative 
uncertainty estimation. It is only necessary to give a 
complete description of the estimation procedure and 
any relevant approximations. 
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